Persistent spin currents in mesoscopic Heisenberg rings 
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We show that at low temperatures T an inhomogeneous radial magnetic field with magnitude B 
gives rise to a persistent magnetization current around a mesoscopic ferromagnetic Heisenberg ring. 
Under optimal conditions this spin current can be as large as 5r/XB(r/ft) exp[— 27r((ji/iBB/A)^''^], as 
obtained from leading-order spin-wave theory. Here g is the gyromagnetic factor, /ib is the Bohr 
magneton, and A is the energy gap between the ground state and the first spin-wave excitation. 
The magnetization current endows the ring with an electric dipole moment. 
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The controlled fabrication of submicron devices has 
opened the door to a rich new field of theoretical and 
experimental physics. At low temperatures these devices 
are mesoscopic in the sense that their quantum states 
must be described by coherent wave functions extending 
over the entire system. Then the usual assumptions un- 
derlying the averaging procedure in statistical mechanics 
are not necessarily valid, and quantum-mechanical inter- 
ference effects become important Q. 

A prominent example is persistent currents in meso- 
scopic normal metal rings threaded by a magnetic flux 

S. Although this phenomenon was predicted long ago 
01 , the experimental difficulties in measuring persis- 
tent currents in an Aharonov-Bohm geometry were only 
overcome in the past decade j3, IS B| ■ Surprisingly, for 
metallic rings in the diffusive regime the observed cur- 
rents were much larger than predicted by theory 
On the other hand, in the ballistic regime |^ the order 
of magnitude of the observed current can be explained 
with a simple model of free fermions moving on a ring 
pierced by a magnetic flux (p. Then the stationary en- 
ergies are e„ = fi^A:^/2m*, where fc„ = T"('^ ~ 
n = 0,±l,±2,..., are the allowed wavevectors for a ring 
with circumference L. Here — hc/e is the flux quan- 
tum and TO* is the effective mass of the electrons. In 
the simplest approximation, one may calculate the cur- 
rent / — —cdVl^c{4>)/d<j) at constant chemical potential 
/i from the flux-dependent part of the grand canonical 
potential fJgc (</>). At finite temperature T, one obtains 
for spinless fermions 



.(e„-A.)/T ^ I ' 



(1) 



where Vn — hkn/m^. For T <ti jj, the amplitude of the 
current is /max ~ —evp/L (where vp\s the Fermi veloc- 
ity), in agreement with experiment [g- 

In this Letter, we show that Heisenberg spin chains 
in inhomogeneous magnetic fields can be used to real- 
ize a spin current analogue of mesoscopic persistent cur- 
rents in normal metal rings. Note that in the presence 
of spin-orbit coupling spin currents in spin chains can 



also be driven by inhomogeneous electric fields due 
to the Aharonov-Casher effect 0. As detailed later on, 
the magnetization current is carried by magnons and en- 
dows the ring with an electric dipole field, which is the 
counterpart of the magnetic dipole field associated with 
the persistent charge current in a normal metal ring. We 
find that for realistic parameters the spin analogues of 
the experiments in Refs. 4, 5,, ,6 require the detection of 
a potential drop on the order of nanovolts. 

Due to its relevance for information processing based 
on spin degrees of freedom, the problem of magnetization 
transport has received a lot of attention recently ^3 > es- 
pecially for systems where spin currents are carried by 
itinerant electrons IH El El El El E3, E3 • Magnetic 
insulators also show interesting spin transport phenom- 
ena Very recently, Meier and Loss [l9| calculated 
the mesoscopic spin conductance for Heisenberg-type sys- 
tems in a two-terminal geometry. Here we consider the 
same problem in a ring geometry and for inhomogeneous 
magnetic fields that lead to noncoplanar spin configura- 
tions. It is known that such configurations can lead to 
dissipationless transport of charge and spin [ill Q, ^1 . 
In the following, we use linear spin- wave theory to derive 
the mesoscopic persistent spin current circulating in the 
ring, corresponding to an infinite spin conductance, and 
explicitly obtain the associated electric field. 

Let us start with a general Heisenberg Hamiltonian 



H ^ - ^ Jij Si ■ Sj 



5MB 



Bi ■ Si 



(2) 



where the sums are over all sites of a chain with peri- 
odic boundary conditions, Jij are general exchange cou- 
plings, and Si are spin-5 operators normalized such that 
Sf = S{S + 1). The last term in Eq. @ is the Zee- 
man energy associated with an inhomogeneous magnetic 
field Bi — B{ri). We assume that the magnetic field at 
each lattice site is sufficiently strong to induce permanent 
magnetic dipole moments = g^^{Si), not necessar- 
ily parallel to Bi, where (. . .) denotes the usual thermal 
average. Moreover, we assume that the unit vectors 
= mi/\mi\ trace out a finite solid angle Q, on the unit 



FIG. 1: Classical spin configuration rhi of a nearest-neighbor 
ferromagnetic Heisenberg ring in a radial magnetic field Bi. 



sphere in order-parameter space as we move once around 
the chain. The simplest geometry is a ferromagnetic ring 
in a crown-shaped magnetic field, as illustrated in Fig.^ 

In the general case, the Hamiltonian ^ implies the 
equation of motion 



OSi \ — ^ 



0, 



(3) 



where hi — gfisBi, and li^j — JijSi x Sj is the spin 
current from site to site rj. From Eq. it is easy 
to show that in equilibrium r hj ■ (li^j) — which is 
Kirchhoff's law for spin currents |20|. 

It is convenient to decompose the spin operators as 
Si = sfrhi + S^, with S^ ■ rhi ~ 0. Substituting this 
into Eq. Q we obtain H = H^^ + + H' , with 



H' 



2 Z^-y 



J a • S A 



h, 



J2hr-ni,S^^ , (4) 

i 

(5) 
(6) 



To develop the standard semiclassical spin-wave theory, 
we expand in the small parameter 1/5'. The leading- 
order term [i.e., Eq. Q with sf replaced by S] yields the 
classical energy. A necessary condition for its minimum 
is the invariance under small variations Srhi. This leads 
to the condition 



h, 



5^ Jijthj 



0, 



(7) 



which shows that the magnetization aligns parallel to the 
sum of external and exchange field, as expected. For 
given hi and Jij, this is a system of non-linear equa- 
tions for the spin directions rhi in the classical ground 
state. To study fluctuations, we expand the transverse 
components as S:^ = ^ J2p=± ^i^^^i^ with spherical ba- 
sis vectors = ej + ipef, p = ±, where {ej,ef,rhi] is 



FIG. 2: Definition of the rotated basis vectors and the corre- 
sponding rotation angles uji^j and uoj^i. The rotated vectors 
el = ej lie in the intersection of the two tangent planes asso- 
ciated with rhi and rhj, shown as a dashed line. The dotted 
line is the geodesic connecting the two sites on the unit sphere 
in order-parameter space. 



a local orthogonal triad of unit vectors. The transverse 
part of our spin Hamiltonian can then be written as 



(8) 



P,! 



tionSf ^S-b\bi and 5+ 



For explicit calculations, we use the standard representa- 

(5r)t ^ V25fe, [1 + 0{S-1)] 
in terms of canonical boson operators hi. Note that 
= 0{S), whereas H' = OiS^'"^) due to Eq. CZJ). 

While the rhi are fixed by Eq. (TJ , there is a remaining 
local U{1) gauge freedom associated with the rotation of 
the transverse basis vectors e^. Let us rewrite the scalar 

product ef • such that the local gauge invariance is 
manifest. Rotating the transverse basis vectors with an 
angle oJi^j around the local normal rhi leads to — 
^ipuji^j^p^ as shown in Fig. |2| If we choose the rotated 
basis vectors and associated with two neighboring 
points rhi and rhj on the unit sphere such that e| and 
are equal and parallel to rhi x hij, then e^-e^ = 
{rhi ■ hij — pp'). Rotating back to the original basis, we 
arrive at 



(rhi ■ fnj — pp') eKp[ipLUi 



- ip LOj 



(9) 



Note that the condition ef = e| implies Im[e^ j 
0. Geometrically, this means that and are related 
by parallel transport along the shortest path (i.e., 
a geodesic, corresponding to a rotation around an axis 
parallel to rhi x hij) connecting the points rhi and rhj 
on the surface of the unit sphere. Substituting Eq. 
into Eq. (jSJ, we obtain 

= ^ E -^^^ [(1 + ■ rh,)e^^-^-^-^-^^SrS+ 

- (1 - • rhj)e'^'^^-^+'^^-'^ S- + h.c.l . (10) 
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This expression is manifestly invariant under the local 
U{1) gauge transformation uji^j — > tUi^j + ai, Sf 
gPgipai^ with Ui arbitrary. The first term in Eq. H1U|) 
shows that a localized spin deviation acquires a phase as 
it moves between sites, which is due to the site-dependent 
orientation of the ground-state magnetization. The local 
gauge invariance implies conservation of the associated 
current, 

3 i 

To see this explicitly, we use the equation of motion 
to write 



-St • [rh, X (/i, - ^ J^3S]rn,^ ] ■ (12) 



we shall assume a radial magnetic field Bi = |B|ri/|ri|, 
with the spins located at constant latitude %% = t? (see 
Fig. HI). For \h\ = g^lB.\B\ > JS{2TT/Nf, the classical 
ground-state configuration rh; is radial as well, with a 
slightly different latitude dm satisfying 

sin(t9„ - i9) = ~{JS/\h\) [1 - cos(27r/iV)] sin(2i3„) . 

^ (14) 

At low temperatures (T ^ JS) we may approximate 
the magnon energies by e„ -I- \h\ with dispersion e„ — 
JSa^k^, where a is the lattice spacing and kn = ^{n ~ 
^) are the quantized wavevectors of the magnons on a 
ring with circumference L. To leading order in spin-wave 
theory, we then obtain the following from Eq. H13|) for 
the magnetization current = {g^^/h)Ig: 



L 



,(£„ + |h,|)/T _ I 



(15) 



Within linear spin-wave theory, the last term may be ne- 
glected since it is an order S'""'^/^ smaller than the first 
term, due to the condition Q. Taking the thermal av- 
erage of both sides of Eq. (|12() and using the fact that 
equilibrium averages arc time- independent, we conclude 
that the average of Eq. indeed vanishes, correspond- 
ing to a longitudinal spin current with vanishing lattice 
divergence. 

However, in a ring geometry there can be a finite cir- 
culating spin current in thermal equilibrium provided 
the classical spin configuration rrii covers a finite to- 
tal solid angle f2 on the unit sphere in order-parameter 
space as we move once around the ring. In the case 
of nearest-neighbor coupling, we obtain the explicitly 
gauge-invariant expression Vl = X^iLi ~ ^i^i-i], 
where N is the number of spins. Note that Vl can be 
identified with the total defect angle ( "anholonomy" ) as- 
sociated with the corresponding parallel transport of a 
tangential vector along a closed path of geodesies . 

We now evaluate the spin current for a ferromag- 
netic Heisenberg chain with nearest-neighbor coupling 
Ji,i±i = —J < 0. The component of the spin current 
in the direction rh; , the divergence of which appears in 
Eq. H12|) . can be written in the gauge invariant form 

Is = - J(m, • {St X St+i)) = -dFsin)/dn , (13) 

where Fs (il) is the free energy of the spin system. On the 
other hand, the transverse spin current component in the 
direction of rhi x rhj is an order 1 /N smaller than Is . In 
deriving Eq. (|13|l . we have neglected the terms involving 
the combinations SJ' and 3^3^ in Eq. (^U)), because 
for a ferromagnet the contribution of these quantum fluc- 
tuations to Is involves higher powers of 1/A^ which are 
dominant only for T — )■ 0. Since the current is a topolog- 
ical property of the system, we may choose any conve- 
nient geometry for explicit calculations. For simplicity. 



where u„ = h ^den/dkn — 2JSa^kn/h. Clearly, Eq. ifT^ 
is the exact bosonic analogue of Eq. |^ • 

If the temperature is large compared with the level 
spacing A — JS{2it/N)'^ between the ground state and 
the first magnon excitation, then Eq. (|15|l can be evalu- 
ated analytically. 



sinfi 



h cosl7-cosh(27rVlWA) 

(16) 

which is accurate for A <c; T ^ JS. In view of 
|/i|/A = S{giiB.\B\/JS'^){N/2TTf, it is clear that the 
persistent spin current is a mesoscopic quantum effect, 
which vanishes for an infinite system or in the classical 
limit S' — > oo with constant JS^. The radial spin con- 
figuration shown in Fig. ^ yields SI = 27r(l — cosdm) to 
leading order in 1/iV, so that the factor | sinfij is maxi- 
mal for cos?9„ = 3/4 or 1/4, i.e., -0^ ~ 41° or 76°. For 
|/i| w A Eq. IHI) predicts a magnetization current on the 
order of Im ~ g^BT/h around the ring. 

Experimentally, the magnetization current can be de- 
tected by measuring the electric voltage between two 
points above and below the ring. The source of a spa- 
tially varying magnetization M{r) = J^i^i''" ~ '"j)'^* 
is an effective current density j{r) = cV x M{r). If 
the magnetization is moving with velocity v{r), it is ac- 
companied by a polarization P{r) = x M{r) to 
leading order in v(r)/c. This is easily shown by means 
of a Lorentz boost to the rest frame of the magnetic 
dipoles ,23j- The polarization corresponds to a charge 
density p{r) = — V • P(r), which in turn generates an 
electric field E{r) = — V0(r), with the scalar potential 
(j>{r) = J d^r' p{r')/\r — r'\. Combining the above re- 
lations, we obtain a generalized Biot-Savart law for the 
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FIG. 3; Lines of constant electric potential due to a spin cur- 
rent circulating in a mesoscopic Heisenberg ring with radius 
R in the geometry of Fig. Qfor i?™ = 30°. 

scalar potential due to magnetic dipole currents, 
^(r) = \ j d^'Hr') X M{r')] ■ 

- 'ffidr'.rhir')]-j^^. (17) 

The second line is valid for a current loop, where 
v\M\d^r — I„idr. Although Eq. ((T7|l is of fundamen- 
tal importance in spin transport, we have not been able 
to find it in standard texts on classical electrodynam- 
ics. The expression derived in Ref. for a straight 
line is a special case of Eq. (|17|) . For the simple ge- 
ometry shown in Fig. ^ the integration in Eq. H17(l can 
be reduced to elliptic integrals. The equipotentials are 
shown in Fig. O In the far zone, the electric field ap- 
proaches a dipole field, with potential (j){r) = p ■ r/|rp 
and dipole moment p = —ez{Im/c)Lsmdm- To esti- 
mate the order of magnitude of the magnetization cur- 
rent, consider a mesoscopic S = 1/2 Heisenberg chain 
with g = 2, N = 100, and J = lOOK. Then the con- 
dition gnB\B\ w A is satisfied for \B\ O.IT. To ob- 
tain a sizable one should generate inhomogeneities of 
the magnetic field in the submicron range; these may be 
achievable in the vicinity of a magnetic flux line trapped 
in a type-II superconductor. In the dipole approxima- 
tion, the potential drop between two points located a dis- 
tance d above and below the loop on the z-axis is given 
hy U ^ G.5nV • (r/K)(L/nm)/(d/nm)2. For T = 60K 
and d = L = lOOnm this yields a voltage U « 0.3nV. 
However, its experimental detection is difficult, because 
mobile charges will tend to screen this static dipole field. 
Note that the magnetization current in Eq. H15(l involves 
only a single factor of g^J-B^ whereas for a two-terminal 
geometry the current is proportional to {gfi^)'^ [l9| . 



In summary, we have shown that a ferromagnetic 
Heisenberg ring in an inhomogeneous magnetic field can 
support a persistent magnetization current, which is the 
precise bosonic analogue of the persistent charge current 
in normal metal rings. The magnetization current is a 
mesoscopic quantum interference effect and flows with- 
out dissipation, corresponding to an infinite spin con- 
ductance. For weak magnetic fields {gfiB\B\ « A) the 
magnetization current gives rise to an electric dipole field. 
For larger fields, the ballistic contribution to the spin cur- 
rent considered in this work is exponentially suppressed. 
By analogy with the persistent charge current in nor- 
mal metal rings, we expect that in this regime the per- 
sistent magnetization current is dominated by collective 
phenomena such as spin diffusion and weak localization 
effects. The calculation of the dominant contribution in 
this regime is still an open problem. With suitable mag- 
netic fields, persistent magnetization currents should also 
exist in antiferromagnetic or ferrimagnetic spin chains, 
where we expect finite currents at T = due to quantum 
fluctuations. 

This work was supported by the DFG via Forscher- 
gruppe FOR 412, Project No. KO 1442/5-1. 



[1] Y. Imry, Introduction to Mesoscopic Physics (Oxford 

University Press, Oxford, 1997). 
[2] F. Hund, Ann. Phys. (Leipzig) 32, 102 (1938). 
[3] M. Biittiker, Y. Imry, and R. Landauer, Phys. Lett. 96A, 

365 (1983). 

[4] L. P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, 

Phys. Rev. Lett. 64, 2074 (1990). 
[5] V. Chandrasekhar, R. A. Webb, M. J. Brady, M. B. 

Ketchen, W. J. Gallagher, and A. Kleinsasser, Phys. Rev. 

Lett. 67, 3578 (1991). 
[6] D. Mailly, C. Chapelier, and A. Benoit, Phys. Rev. Lett. 

70, 2020 (1993). 
[7] H.-F. Cheung, Y. Gefen, and E. K. Riedel, IBM J. Res. 

Develop. 32, 359 (1988). 
[8] Z. Cao, X. Yu, and R. Han, Phys. Rev. B 56, 5077 (1997). 
[9] Y. Aharonov and A. Casher, Phys. Rev. Lett. 53, 319 

(1984). 

[10] See, for example, D. D. Awschalom, D. Loss, and N. 
Samarth (Eds.), Semiconductor Spintronics and Quan- 
tum Computation, (Springer, Berlin, 2002). 

[11] D. Loss, P. Goldbart, and A. V. Balatsky, Phys. Rev. 
Lett. 65, 1655 (1990). 

[12] A. Stern, Phys. Rev. Lett. 68, 1022 (1992). 

[13] X.-C. Gao and T.-Z. Qian, Phys. Rev. B 47, 7128 (1993). 

[14] J. Konig, M. C. Bonsager, and A. H. MacDonald, Phys. 
Rev. Lett. 87, 187202 (2001). 

[15] G. Tatara and H. Kohno, Phys. Rev. B 67, 113316 
(2003). 

[16] A. G. Mal'shukov, C. S. Tang, C. S. Chu, and K. A. 

Chao, cond-mat/0211559 
[17] S.-Q. Shen and X. C. Xie, Phys. Rev. B 67, 144423 

(2003). 

[18] L. Y. Gorelik, R. I. Shekhter, V. Vinokur, D. Feldman, 



5 



V. Kozub, and M. Jonson, cond-mat/0211563 
[19] F. Meier and D. Loss, Phys. Rev. Lett. 90, 167204 (2003). 
[20] Note that the spin current is a second rank tensor. We call 

riii-Ii^j the longitudinal spin current, where longitudinal 

refers to spin space. 
[21] See, for example, A. Shapere and F. Wilczek, Geometric 

Phases in Physics, (World Scientific, Singapore, 1989). 



[22] To compare Eqs. Q and (I15II . we should divide the 
currents by the corresponding "charges", i.e., the par- 
ticle number current //(— e) should be compared with 

Is/h = Im/igfj-s)- 

[23] J. E. Hirsch, Phys. Rev. B 60, 14787 (1999). 



